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Abstract 

Thermospheric infrared radiance at 4.3 /jrn is susceptible to the influence of solar- 
geomagnetic disturbances. Ionization processes followed by ion-neutral chemical re- 
actions lead to vibrationally excited NO + (i.e. , NO + (v)) and subsequent 4.3 /irn 
emission in the ionospheric E-region. Large enhancements of nighttime 4.3 gm emis- 
sion were observed by the TIMED/SABER instrument during the April 2002 and 
October-November 2003 solar storms. Global measurements of infrared 4.3 gm emis- 
sion provide an excellent proxy to observe the nighttime E-region response to au- 
roral dosing and to conduct a detailed study of E-region ion-neutral chemistry and 
energy transfer mechanisms. Furthermore, we find that photoionization processes 
followed by ion-neutral reactions during quiescent, daytime conditions increase the 
NO + concentration enough to introduce biases in the TIMED/SABER operational 
processing of kinetic temperature and CO 2 data, with the largest effect at summer 
solstice. In this paper, we discuss solar storm enhancements of 4.3 gm emission ob- 
served from SABER and assess the impact of NO + (v) 4.3 gm emission on quiescent, 
daytime retrievals of Tk/CC >2 from the SABER instrument. 
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1 Introduction 


The Thermosphere-Ionosphere- Mesosphere Energetics and Dynamics (TIMED) 
satellite was launched in December 2001. The primary objective of TIMED 
is to investigate and understand the energetics of the mesosphere, lower ther- 
mosphere and ionosphere (MLTI) region. Precipitating energetic particles in- 
duced by solar-geomagnetic disturbances rapidly and dramatically change the 
chemistry and thermal energy balance of the MLTI region. Several major 
solar eruptive events have occurred during TIMED mission operations. The 
two events considered in this paper are the April 2002 and the Halloween 
(October-November) 2003 solar storms. TIMED provides a suite of observa- 
tions uniquely suited for investigating the response of the MLTI region to 
solar-geomagnetic disturbances. For example, data from the TIMED Global 
Ultraviolet Imager (GUVI) instrument have been used to analyze the total en- 
ergy flux and characteristic energy of auroral electrons during the April 2002 
storm (Christensen et aL, 2003) and to investigate the causal mechanisms 
responsible for nighttime detached auroras during the Halloween 2003 storm 
(Zhang et aL, 2005). 

The Sounding of the Atmosphere using Broadband Emission Radiometry 
(SABER) instrument is an infrared limb sounder on the TIMED satellite 
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(Russell et al., 1999). SABER provides global measurements of broadband 
infrared limb emission, which are analyzed to derive vertical profiles of kinetic 
temperature and key chemical species needed to quantify the radiative and 
chemical energy balance of the MLT1 region. Large enhancements of thermo- 
spheric infrared emission were observed by several of the SABER radiometer 
channels during the April 2002 and Halloween 2003 solar storm events. Radi- 
ance enhancements observed in the 5.3 /mi channel are due to emission from 
the vibration-rotation bands of nitric oxide (NO). The NO 5.3 /irn emission is 
indicative of the conversion of solar energy to infrared radiation and represents 
a “natural thermostat” by which heat and energy are efficiently lost from the 
thermosphere to space and to the lower atmosphere. SABER 5.3 /irn emission 
observations during the April 2002 storm and their subsequent interpretation 
in terms of energy loss were reported by Mlynczak et al. (2003, 2005). 

Another thermospheric infrared spectral region subject to solar-geomagnetic 
influences is the 4.3 /irn region, due to the vibration- rotation bands of NO + 
(Picard et al., 1987; Espy et al., 1988). During solar-geomagnetic storms, elec- 
tron precipitation increases the ionization of the neutral atmosphere produc- 
ing primarily N^~, O^, 0 + and N + (e.g., Banks et ah, 1974; Strickland et al., 
1976). In the ionospheric E-region, these ions react with neutral species to 
produce NO + (Torr et ah, 1990; Fox and Sung, 2001). Some of the ion-neutral 
reactions are exothermic enough to produce vibrationally excited NO + , i.e., 
NO + (v) (Winick et al., 1987). The exothermic reactions are relatively fast. 
Above ~ 110 km quenching of NO + (v) becomes less important and prompt 
emission of NO + (v) at 4.3 /irn is a direct indication of auroral dosing in the E- 
region. The NO + (v) fundamental band (2344 cm -1 ) is nearly coincident with 
the strong major isotopic CC^rq) band at 2349 cm -1 . Infrared emission at 4.3 
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ji m during strong nighttime aurora can be enhanced by several orders of mag- 
nitude by NO + (v) compared to the background CC^rq) radiative emission. 
The E-region is largely inaccessible to in-situ observations and the SABER 
4.3 /mi measurements provide an excellent dataset to investigate the E-region 
response to solar-geomagnetic storms. 

Infrared NO + (v) emission has important implications for the operational pro- 
cessing of SABER routine data products at high-latitudes near summer sol- 
stice, even for quiescent conditions. During daytime, kinetic temperature (Tk) 
and carbon dioxide (CO 2 ) volume mixing ratio (vmr) are simultaneously re- 
trieved in version 1.06 data processing from SABER measurements of 15 /mi 
and 4.3 /im limb emission, respectively (Mertens et ah, 2002). Nighttime Tk 
is retrieved from 15 /mi measurements using CO 2 data from the TIME-GCM 
climatology (Mertens et ah, 2001, 2004). Solar EUV photons ionize the neutral 
atmosphere during the daytime, eventually leading to NO + (v) and emission at 
4.3 /rm through the exothermic ion-neutral reactions mentioned above. These 
processes are not included in the SABER operational Tk/CCD retrieval algo- 
rithm. Quiescent daytime thermospheric NO + (v) 4.3 /mi emission is a small 
percentage of the background 002 (^ 3 ) emission. However, 4.3 /mi radiation 
transfer is strongly non-linear along the limb line-of-sight at mesospheric tan- 
gent heights, with non-negligible contributions coming from the lower thermo- 
sphere. Consequently, weak atmospheric emission from NO + (v) not accounted 
for in the Tk/CC >2 retrieval algorithm can have a rather significant effect on 
retrieved C0 2 throughout the MLTI region, which subsequently introduces a 
bias in retrieved Tk. Most of the SABER data products require temperature as 
input into their respective retrieval algorithms. Thus, a bias in Tk introduces 
biases in nearly all of the SABER data products. We have evidence, from 
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a combination the SABER CO 2 data and model simulations, that strongly 
suggests that quiescent 4.3 pm emission by NO + (v) is introducing errors in 
retrieved Tk/C 02 at high latitudes, with the largest effect at summer solstice. 

The focus of this paper is on NO + (v) 4.3 urn emission during both quiescent 
and geomagnetically disturbed conditions. We present SABER observations of 
nighttime auroral enhancements of 4.3 pm emission during the April 2002 and 
Halloween 2003 storms. We conduct a preliminary assessment of the impact of 
quiescent daytime NO + (v) 4.3 pm emission on high-latitude SABER Tk/CCH 
retrievals. 


2 Solar Storm Enhancements of 4.3 pm Emission 

Nighttime enhancements of 4.3 pm emission can be several orders of magnitude 
during strong aurora. Figure 1 shows representative radiance profiles measured 
by the SABER (channel 7) 4.3 pm radiometer channel during the April 2002 
and Halloween 2003 solar storms. The vertical, dashed line in the figure is the 
channel 7 noise equivalent radiance (NER), which is 7.35 x 10 - ' W/m 2 /sr. 
Prior to the onset of the storm periods, the radiance decreases with increasing 
altitude reaching the detector noise level by ~ 135 km. During the peak of the 
solar storms, the 4.3 pm radiance is enhanced by up to several orders of mag- 
nitude and doesn’t reach the detector noise level until well above 180 km. The 
storm-enhanced radiance profile during the Halloween 2003 event in Figure 1 
has the same slope as its corresponding pre-storm radiance profile, while this 
is not the case for the April 2002 storm event. The correlation between the 
slope of the storm and pre-storm radiance profiles during the Halloween event 
is fortuitous, since the vertical profile of the storm-enhanced radiance profile 
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depends on the specific mapping of the auroral electron energy spectrum on 
the ionization rate profile and subsequent ion-neutral compositions. 

Figure 2 shows a comparison of nighttime precipitating total electron energy 
flux with SABER-derived NO + (v) Volume Emission Rate (VER) at 110 km 
in the northern hemisphere. The NO + (v) VER is the quantity that charac- 
terizes the storm-time enhancement of 4.3 /jrn emission (see section 3.1 and 
Mertens et al. (2007a, 2007b)). The total electron energy flux is determined 
from measurements made by the National Oceanic and Atmospheric Adminis- 
tration (NOAA) Total Energy Detector (TED) instruments onboard the Polar 
Orbiting Environmental Satellites (POES). The NOAA/POES satellites are 
in sun-synchronous orbits. POES-17 is in an orbit that crosses the equator 
northbound at 2215 local time at the sub-satellite point. As a result, POES- 
17 covers the northern hemisphere auroral zone in the time sector ranging 
from about 2000 local time through 0100 local times. POES-15 and POES- 
16 cross the equator at about 1900 and 1400 local times at the sub-satellite 
points. Thus, they sample the northern hemisphere largely during the day- 
time hours. The satellite coverage described above is opposite in the southern 
hemisphere. The total electron energy flux is derived from POES directional 
energy flux measurements over the energy range from 50 eV to 20 keV. Au- 
roral electrons with these energies deposit their energy at altitudes above 100 
km (Rees, 1989). 

The SABER and POES data in Figure 2 were taken during the Halloween 
2003 storm event on 30 October. The data are presented in polar magnetic 
coordinates (magnetic latitude and magnetic local time (MLT)). A contin- 
uous distribution of total electron energy flux was obtained by binning and 
averaging the POES data on a 1x2x3 magnetic latitude (degrees), magnetic 
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longitude (degrees), and UT-time (hours) grid. The gridded POES data are 
interpolated to fill in all magnetic longitude grid points. The corresponding 
MLT are determined by the geocentric-magnetic coordinate transformation. 
The SABER NO + (v) VER are shown at 110 km for each orbit that falls within 
the 3-hour UT-time intervals shown for the corresponding POES data. Au- 
roral precipitation changes rapidly in space and time. Nevertheless, one can 
observe from Figure 2 that high (low) levels of auroral precipitation generally 
correspond to high (low) levels of NO + (v) VER. 


Thus, thermospheric infrared emission at 4.3 /irn is an excellent proxy for 
observing the ionospheric E-region response to strong nighttime auroral dos- 
ing. The combination of SABER 4.3 /irn radiance measurements, observations 
of precipitating electron energy characteristics, ionospheric chemistry and ra- 
diation transfer models, can be used to conduct a detailed investigation of 
our understanding of the physical mechanisms responsible for the E-region 
energetics and the response to auroral electron dosing, since many of these 
mechanisms and their subsequent rates of reactions are still largely uncertain 
(Winick et al., 1987). This work is currently under investigation (Mertens et 
ah, 2007a, 2007b, 2008b). In the next section we briefly summarize our NO + (v) 
energetics model. In section 4 our model simulations are used to assess the 
impact of quiescent, daytime NO + (v) 4.3 /irn emission on the SABER Tk/CCU 
version 1.06 retrievals for high-latitude summer solstice conditions. 



3 NO + (v) Energetics Model 


3.1 NO + (v) Kinetics and Chemistry 


A number of ion-neutral reactions are sufficiently exothermic to produce vi- 
brationally excited NO + (Winick et ah, 1987; Caledonia et ah, 1995; Dothe 
et ah, 1996; Smith et ah, 2000). The reactions that have received the most 
attention are as follows: 


N+ + 0 
N + + 0 2 


o + + n 2 

Q+ + NO 


NO + (u) 

NO + (u) 

NO + (u) 

NO + (u) 

NO + (u) 

NO + (u) 

NO + (u) 

NO + (u) 

NO + (u) 


+ N( 4 S) + 3.08 eV (i'max = 10) 
+ N ( 2 D) + 0.69 eV (i'max — 2) 
+ 0( 1 S) + 2.43 eV (i'max = 8) 

+ 0( 4 D) + 4.66 eV ( i'max — 18) 
+ 0( 3 P) + 6.63 eV (i'max = 28) 
+ N( 4 S) + 1.1 eV (i'max = 3) 

+ 0 2 (X) + 2.83eV (i'max — 9) 

+ 0 2 (a ) + 1.85 eV (i'max = 6) 

+ 0 2 (b ) + 1.20 eV (i'max = 4). 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


Loss of NO + (v) is by collisional quenching and spontaneous radiative emission. 
The quenching and radiative loss processes are represented by 


NO + (u) + N 2 (v = 0) ^ NO + (u - 1) + N 2 (v = 0) 
NO + (v) -»• NO + (v - 1) + hv. 


( 5 ) 

( 6 ) 


Reaction (2) is the only one for which the vibrational distributions of the 
products have been measured in the laboratory (Smith et ah, 1983; denoted 
SM83 in Table 1). Albritton et ah (1979) and Langford et ah (1986) determined 
the branching ratios for electronically excited atomic oxygen atoms formed 
by reaction (2), and Smith et ah deduced the fractional rate that populates 
v — 1 — 14. Reaction (2) is an important contributor to NO + (v) 4.3 /mi 
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emission below ~ 140 km (Winick et al., 1987; Caledonia et al., 1995). In our 
analysis, we assume the total reaction rate and fractional vibrational state 
distribution are those obtained by the above measurements. 

Above 140 km, reaction (1) dominates the production of NO + (v) and subse- 
quent 4.3 /jin emission (Winick et al., 1987; Caledonia et al., 1995; Dothe et 
al., 1996; Smith et al., 2000). However, the vibrational distribution of NO + (v) 
from (1) has never been measured in the laboratory. On the other hand, mod- 
eling of NO chemistry and analysis of high- resolution atmospheric emission 
spectra suggest that N( 2 D) is the dominant channel in reaction (1) (Winick 
et al., 1987; Dothe et al., 1996). There are a number of assumptions for the 
vibrational state distribution reported in the literature. We assume the nom- 
inal case utilized by Winick et al. (1987) (denoted WK87 in Table 1), where 
the vibrational state distribution is given by f(/J=0) = 0.25, f(u=l) = 0.50, 
and f(u=2) = 0.25. For the total rate of reaction we use the rate coefficient 
measured by McFarland et al. (1974) (denoted MF74 in Table 1). 

Reaction (3) has received attention recently as a candidate for non-thermal 
rotational excitation of NO + (z/,J) (Smith et al., 2000; and Duff and Smith, 
2000). Smith et al. (2000) analyzed CIRRIS 1A data and showed that most 
of the NO + produced by reaction (3) is rotationally thermalized below 150 
km. Furthermore, reaction (3) is not a dominant contributor to the overall 
production of NO + (v) below ~ 160 km. Model studies conducted by Duff and 
Smith (2000) showed that NO + (v) atmospheric radiance from reaction (3) 
was dominated by v =1 emission. Taking these results into consideration, and 
restricting our analysis to altitudes below ~ 160 km, we assume rotational LTE 
for this reaction and assume that NO + is created in vibrational state v = 1 
(denoted DS00 in Table 1). We also use the reaction rate coefficient measured 
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by McFarland et al. (1973) (denoted MF73 in Table 1), as suggested by Smith 
et al. (2000). 

The large exothermicity of reaction (4) makes the production of NO + (v) possi- 
ble. This process is important for high [NO], which occurred during the recent 
solar storm events (Richards, 2004; Mlynczak et al, 2003). However, there were 
no previous observation-based studies that presented evidence of NO + (v) ex- 
citation and subsequent 4.3 /irn emission from this reaction. On the contrary, 
in a recent study we found this reaction to be the dominant source of NO + (v) 
excitation and auroral 4.3 /jrn emission below 160 km (Mertens et ah, 2008b). 
The total rate coefficient for reaction (4) is obtained from the measurements 
of Lindinger et al. (1974) (denoted LD74 in Table 1). The nascent NO + (v) 
distribution from this reactions is given by Mertens et al. (2008b) (denoted 
MT08 in Table 1). However, since we are considering quiescent conditions in 
this paper, we ignore this reaction in the subsequent analysis. 

Collisional quenching of NO + (v) in (5) is negligible above 130 km, but becomes 
an important loss mechanism by 110 km (Winick et ah, 1987). Doblcr et 
ah (1983) (denoted D083 in Table 1) and Federer et ah (1985) measured 
similar rate coefficients for the quenching of NO + (^=l) in neutral collisions. 
Quenching rates of levels with v > 1 have not been measured. As a result, 
various quenching schemes have been assumed in the literature. We use the 
scheme of Winick et ah (1987), where quenching of NO + (v) is by collisions 
with N 2 and single-quantum quenching is assumed for all NO + (v) levels. 

Radiative loss of NO + (v) in (6) is modeled using the Einstein A-coefficients 
reported by Werner and Rosmus (1982), who calculated spontaneous emission 
rates for the fundamental and first two overtone emission sequences for v = 
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1 to v = 20. We use the fundamental (Au—1) and the first overtone (Az/=2) 
A-coefficients, as the second overtone spontaneous emission rates are much 
smaller than the Au=l and Av=‘l rates. 

The reactions rate coefficients, vibrational state distributions, collisional and 
radiative loss rates described in this section allow one to calculate the steady- 
state NO + vibrational state populations, provided one knows the ion and 
neutral concentrations that appear on the left-hand side of reactions (l)-(4). 
In our NO + (v) energetics model, the ion/neutral concentrations are calculated 
using the FLIP model, as described in the next subsection. 


3.2 FLIP Model 


The field-line interhemispheric plasma (FLIP) model is a one-dimensional 
ionosphere and thermosphere model that calculates plasma densities and tem- 
peratures along entire magnetic flux tubes from below 100 km in one hemi- 
sphere, through the plasmasphere to below 100 km in the conjugate hemi- 
sphere (Richards, 2002). The equations solved along the flux tubes are the 
time-dependent continuity and momentum equations for 0 + , H + , He + , and 
N + . The electron and ion temperatures are obtained by solving the energy 
equation. Chemical equilibrium densities are calculated for NO + , 0 2 , N 2 , 
0 + ( 2 P), and 0 + ( 2 D) ions below 500 km in each hemisphere. The densities of 
minor neutral species NO, 0( 1 D), N( 2 D), and N( 4 S) are obtained by solving 
the time-dependent continuity and momentum equations from 100 to 500 km 
in each hemisphere. The model also solves for the first five vibrational levels 
of molecular nitrogen (N 2 (v)). 
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We use the FLIP model in the simulations of NO + (v) 4.3 /mi emission. FLIP 
provides the ion/neutral species for the reactants in (l)-(3). Figure 3 shows 
results from FLIP simulations. The ion concentrations needed to calculate 
the NO + (v) populations are shown in the left-hand column. The right-hand 
column are the neutral concentrations. The first row is a simulation for a quiet 
day (8 April) prior to the onset of the April 2002 solar storm. The second row 
is a simulation of the ion/neutral concentrations for a representative day (20 
April) during the peak thermospheric infrared response to the April 2002 
storm, with auroral dosing similar to the total electron energy flux observed 
by the NOAA/POES instruments during this storm. Notice that the NO + 
concentration for the auroral dosing is approximately two orders of magnitude 
greater than the concentration during the undisturbed day. The enhancement 
in the NO + concentration for the simulation in Figure 3 is similar to the 
enhancement of thermospheric 4.3 /mi radiance observed by SABER during 
the solar storms (see Fig. 1). 

Figure 4 shows a comparison of FLIP electron density (i.e., [e]) simulations 
at 115 km with ground-based VHF incoherent scatter radar measurements 
of electron density taken by the European Incoherent SCATer (EISCAT) fa- 
cility at Trornsp, Norway (70N, 19E). We compare FLIP/EISCAT [e] during 
a period of low geomagnetic activity and continuous EISCAT measurements 
from June 14-19, 2004. FLIP [e] simulations are in excellent agreement with 
the 1-hour averaged EISCAT [e] measurements during the daytime hours. The 
largest differences occur over several hours near midnight MLT due to ‘quiet- 
time’ auroral activity. A geomagnetic substorm appears to have occurred on 
June 16 (Akasofu, 1964), which can be associated with particle precipita- 
tion (parameterized by the Hemispheric Power (HP) index), enhanced electric 


13 



fields with both field-aligned and electrojet currents (parameterized by the 
AE- index). However, we are analyzing quiescent, daytime conditions in this 
paper and Figure 4 demonstrates the efficacy of using FLIP to simulate the 
E-region plasma under these conditions at high-latitude. 


3.3 Infrared Radiance Model 

The SABER 4.3 /mi limb emission can be modeled once the NO + (v) pop- 
ulations are computed, which are calculated using the FLIP model and the 
NO + (v) kinetics and chemistry model described in sections 3.2 and 3.1, re- 
spectively. The 002 (^ 3 ) contribution to the SABER 4.3 /irn measurements 
is modeled using the 4.3 /mi channel forward model component of the op- 
erational non-LTE Tk/CCD retrieval algorithm (Mertens et al., 2002). The 
forward model is comprised of two parts: (1) the CO 2 vibrational temperature 
(Tv) model and (2) the limb radiance model. Limb radiance is calculated using 
BANDPAK (Marshall et ah, 1994), which is based on emissivity databases cal- 
culated line-by-line using LINEPAK routines (Gordley et ah, 1994). The CO 2 
Tv model is based on the Modified Curtis Matrix approach (Lopez-Puertas et 
ah, 1986a-b, 1998), and uses BANDPAK in all the radiation transfer calcula- 
tions. 

The in-band SABER NO + (v) 4.3 /irn limb radiance will be simulated using 
LINEPAK, radiative line parameters from the HITRAN 2000 molecular spec- 
troscopic database (Rothman et ah, 2003), the SABER 4.3 /mi channel relative 
spectral response filter, and NO + (v) Tv’s (or equivalently, the NO + (v) popu- 
lations) calculated using FLIP and the NO + (v) kinetics and chemistry model. 
The HITRAN NO + (v) line parameters are derived from the spectroscopic con- 
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stants and parameters calculated by Billingsley (1973), Huber and Herzberg 
(1979), and Werner and Rosmus (1982). Figure 5 is a simulation of infrared 4.3 
/irn limb emission for quiet and disturbed conditions. The disturbed condition 
is for auroral dosing representative of the total electron energy flux measured 
by the NOAA/POES instruments during the April 2002 solar storm. Clearly, 
NO + (v) dominates the 4.3 /irn spectrum. The black curve is the SABER spec- 
tral response filter. The SABER in-band contributions from NO + (v) is 60 to 
80% of the SABER 4.3 /irn channel radiance for the simulation in Figure 5. 
The SABER in-band NO + (v) 4.3 /irn emission is comprised of NO + (Az/=l) 
transition from v — 1 to 4. 

3-4 Auroral Excitation of N 2 (v) 

Auroral electrons vibrationally excite molecular nitrogen by (1) inelastic col- 
lisions between low energy secondary electrons and N 2 (Newton et al., 1974), 
and (2) chemical reactions (Richards, 2002). Reaction (3) is highly depen- 
dent on the vibrational level of N 2 and can significantly increase the effective 
reaction rate in the F2 region. These mechanisms are included in the FLIP 
model. 

Auroral excitation of N 2 (z/) also affects local 4.3 /irn emission from C0 2 (r , 3) 
below about 110 km through V-V exchange between C0 2 (z/ 3 ) and N 2 (z/): 

C0 2 (u 3 ) + N 2 {y) <— ■> C0 2 (v 3 — 1) + N 2 (y + 1) + 18 cm 1 (7) 

Reaction (7) was measured by Inoue and Tsuchiya (1975) for N 2 (zi=l). 
Above 110 km, N 2 (i/) is quenched by collisions with atomic oxygen and locally 
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decouples from CCH^). However, aurora can potentially influence CCH^) 
4.3 fin 1 emission above 110 km through non-local processes due to V-V ex- 
change between 002 (^ 3 ) and ^(z') below 110 km followed by radiative trans- 
port to higher altitudes. 

The influence of auroral electron excitation of N 2 (zz) on 4.3 fim emission is not 
included in the simulations presented in this paper. Our studies have shown 
that auroral CC^^) 4.3 fin 1 emission to be insignificant in comparison to 
auroral NO + (v) 4.3 fin 1 emission, consistent with the results found by O’Neil 
et al. (2007). This is especially true for the quiescent daytime conditions con- 
sidered in this paper. 


4 Influence of NO + (v) 4.3 fim Emission on SABER Tk/C0 2 Re- 
trievals 

Representative SABER version 1.06 monthly zonal mean CO 2 vmr data are 
shown in Figure 6 for winter, equinox, and summer seasons. Comparisons of 
SABER CO 2 vmr with TIME-GCM (Roble, 1995) model-generated clima- 
tological data show that SABER CO 2 vmr is systematically lower than the 
model simulations at all altitudes above roughly 70-75 km (Mertens et ah, 
2008a). This result is consistent with previous CO 2 vmr measurements de- 
termined from atmospheric infrared remote sensing observations (Kaufmann 
et al, 2002; Zaragoza et ah, 2000; Lopez-Puertas et ah, 2000; and references 
therein). However, the puzzling feature in Figure 5 is the meridional gradient 
in the CO 2 concentration near the polar regions not seen in the model data, 
which dramatically increases from equinox to summer season. In fact, there is 
a step-function feature in the July SABER data coincident with the latitude 
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region of the auroral oval. Simulations from our NO + (v) energetics model sug- 
gest that the step-function feature in the C0 2 data is caused by NO + (v) 4.3 
/irn emission, which is not accounted for in the operational SABER non-LTE 
Tk/C0 2 retrieval algorithm. 

The omission of the NO + (v) contribution to the 4.3 /irn radiance in the 
Tk/C0 2 retrieval algorithm underestimates the radiation transfer modeling of 
the SABER 4.3 //rri channel measurement, which in turn results in an overesti- 
mation of C0 2 . Taking the northern hemisphere in Figure 6 as an example, the 
polar region is exposed to increasing periods and intensity of solar illumination 
as the seasons proceed from winter (January) to summer (July). Subsequently, 
the polar region NO + concentrations increase from winter to summer as the 
rate of photoionization and photoelectron concentrations increase (Mertens et 
al., 2008a). The overestimated C0 2 concentration near summer solstice (i.e., 
the July SABER data) is coincident with the auroral oval latitude band due 
to the geometric effects of Earth’s magnetic field on the photoelectrons that 
participate in the E-region chemistry leading to NO + . 

Figure 7 is a simulation of the SABER 4.3 /irn channel measurements with and 
without NO + (v) 4.3 /.irn emission included. The simulated observations are for 
five SABER measurement profiles near 69N on 4 July, 2002. The NO + (v) 4.3 
/rm emission was calculated using the NO + (v) energetics model described in 
section 3. As shown in Figure 7, NO + (v) contributes ~ 1% of the total limb 
radiance near 110 km and between 6-8% at 130 km. These simulations are 
consistent with the NO + (v) contribution to 4.3 /rm emission determined from 
the spectral radiance observations made by the MIPAS (Michelson Interferom- 
eter for Passive Atmospheric Sounding) instrument on the ENVISAT-1 satel- 
lite (Lopez-Puertas, private communications, 2004). This consistency gives us 
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some confidence in our NO + (v) energetics model. 


The NO + (v) contribution to the daytime 4.3 /jrn emission has a significant 
effect on the CO 2 retrieval. Because of the severe non-linearity in the 4.3 /jrn 
radiation transfer, the weighting functions in the mesosphere extend up into 
the lower thermosphere, with non-negligible contributions arising from the 130 
km region (Mertens et ah, 2008a). Thus, the CO 2 retrieval can be overesti- 
mated throughout the mesosphere as a result of relatively weak NO + (v) 4.3 
H m emission contributions above ~ 110 km. 


An overestimation of CO 2 will cause an underestimation in Tk. The low bias 
in Tk was estimated by forcing the non-LTE Tk/C02 retrieval algorithm to 
use the SABER-C0 2 mean profile at 41N. The TIME-GCM model results 
suggest that the CO 2 vmr in the northern polar region for summer solstice 
should be roughly the same as the CO 2 vmr at mid-latitudes. CO 2 validation 
studies (Mertens et al. 2008a) give us confidence in the SABER C0 2 data 
at mid-latitudes. Temperature retrievals in the polar summer region with the 
anomalously large CO 2 profiles (i.e. , Figure 6) versus using the SABER-CO 2 
mean profile at 41 N should give a reasonable assessment of the Tk bias due 
to NO + (v) 4.3 fin 1 emission. The result of this sensitivity study is shown in 
Figure 8. 


The sensitivity study in Figure 8 indicates a low bias in polar summer Tk of 
about 5 K in the upper mesosphere, which is centered in altitude at 80 km. 
Above 100 km, the bias is much larger. The bias in Tk will introduce errors 
in most of the other SABER data products. 
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5 Summary 


Large enhancements of thermospheric infrared emission was observed by SABER 
during the recent superstorm events: the April 2002 and the Halloween 2003 
solar storms. Radiance enhancements in the SABER 4.3 /mi channel are due 
to vibrationally excited NO + , which is an excellent proxy for observing the 
ionospheric E-region response to auroral electron dosing induced by solar- 
geomagnetic disturbances. An NO + (v) energetics model was developed to 
simulate the exothermic ion- neutral reactions that produce NO + (v) and the 
subsequent emission at 4.3 /mi. Efforts are underway to use this model along 
with SABER 4.3 /irn limb radiance observations and electron energy char- 
acteristics measured by the NOAA/POES satellites to investigate E-region 
ion-neutral chemistry and energy transfer mechanisms (Mertens et al., 2007a, 
2007b, 2008b). 

NO + (v) 4.3 /.irn emission during quiescent, daytime conditions introduce biases 
in version 1.06 SABER.-retrieved Tk/CC >2 at high-latitudes in the mesosphere 
and lower thermosphere. The impact of NO + (v) 4.3 /in r emission on retrieved 
CO 2 is maximum for summer solstice conditions, moderate for equinox, and 
negligible for winter conditions. The NO + (v) contribution to SABER 4.3 /mi 
channel measurements for polar summer solstice is 1% at 110 km and 6-8% at 
130 km, consistent with analysis of MIPAS 4.3 /mi spectral radiance measure- 
ments. The NO + (v) contribution to the SABER 4.3 /mi observations results 
in an overestimation of retrieval CO 2 which subsequently induces a low bias 
in retrieved Tk by 5 K near 80 km, with much larger biases above 100 km. 
The temperature biases will introduce errors in most of the other SABER data 
products. 
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Table 1 

NO + (v) Chemistry Model. Reference abbreviations are defined in section 3.1 


Reaction 

Reaction 

Rate Coefficient 

Reference 

Nascent NO + (' 

Number 

Process 

(cm 3 /s) 

(Rate Coeff) 

Distribution 

1 

N++0 

1.4xl0 -1 °(Tj/300) -0 ' 44 Tj < 1500 K 

MF74 

WK87 



5.2xl0 _11 (Tj/300)°' 2 ° Ti > 1500 K 



2 

n++o 2 

2.7xlO -10 

SM83 

SM83 

3 

o++n 2 

1.2xl0 _12 (Tj/300) _1 T, : < 750 K 

MF73 

DS00 



8.0xl0 _14 (Tj/300) 2 Tj > 750 K 



4 

O^+NO 

4.5xlO -10 

LD74 

MT08 

5 

NO+(v)+N 2 

7.0xl0 -12 

D083 

N/A 
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Fig. 1. SABER 4.3 //in limb radiance measurements during the April 2002 and 
Halloween 2003 solar storms. The blue lines are radiances profiles prior to the onset 
of the storms while the red lines are radiance profiles observed during the peaks of 
the storms. The dashed lines are the noise equivalent radiance (NER) level for the 
SABER 4.3 //nr channel. 

Fig. 2. NOAA/POES in-situ measurements of total energy flux (ergs/cnr 2 /s) and 
SABER-derived NO + (v) Volume Emission Rates (VER) (x 10 -9 ergs/cnr 3 /s) in 
the northern hemisphere on October 30, 2003. The data are shown in magnetic 
coordinates (magnetic latitude and magnetic local time (MLT)). The NOAA/POES 
data are gridded as described in section 2 and shown in 3-hour UT-tinre intervals. 
The SABER-derived NO + (v) VER are shown at 110 km for orbits that fall within 
the corresponding 3-hour UT-tinre intervals of the POES data. 

Fig. 3. FLIP model simulation of ion/neutral densities for April 2002, a represen- 
tative quiet day and a day of peak solar-geonragnetic activity. The auroral dos- 
ing parameters used for the simulation shown in the second row are consistent 
with the maximum daily-averaged auroral electron dosing parameters observed by 
NOAA/POES on 20 April, 2002: total electron energy flux is 9 ergs/cnr 2 /s and the 
characteristic energy is 6 keV. 

Fig. 4. Comparison of FLIP electron density (i.e. , [e]) simulations with ground-based 
VHF incoherent scatter radar electron density measurements taken by the EISCAT 
facility at Trornsp, Norway from June 14-19, 2004. The first, second, and third panels 
correspond to the HP-index, AE-index, and kp-index, respectively. The bottom 
panel shows the actual FLIP/EISCAT [e] (cm -3 ) comparisons at 115 km. The grey 
symbols in the bottom panel represent the individual EISCAT measurements (~ 
1-minute average). The blue line represents the EISCAT 1-hour average [e]. The 
red line is the FLIP [e] simulation written out in 1-hour intervals. 



Fig. 5. Simulations of 4.3 /jrn limb emission during quiescent and auroral conditions. 
The auroral conditions are the same as described in figure 3. 


Fig. 6. SABER version 1.06 zonal mean C02(vmr) for July 2002, January and 
September 2003. 


Fig. 7. Radiance simulation of SABER 4.3 fim channel, with and without NO + (v) 
contribution. 


Fig. 8. Temperature retrieval with anomalous polar summer CO 2 profile versus 
temperature retrieval using SABER-CO 2 mean profile at 41N. 
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SABER 4.3um Nighttime Radiance Measurements 

April 2002 October 2003 
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FLIP/EISCAT Plasma Comparisons 
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4.3 um Band Emission: Upper Panel = Quiescent Night: Lower Panel = Night Aurora 
Tangent Height = 140 km Spectral Resolution = 0.05 cm"' 
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SABER 4.3 um Channel Radiance: July 4, 2002) 
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SABER-Derived Mean Temperature Profiles 


